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Abstract 



We investigate the properties of hysteresis cycles produced by a one-dimensional, pe- 
riodically forced Langevin equation. We show that depending on amplitude and fre- 
quency of the forcing and on noise intensity, there are three qualitatively different 
types of hysteresis cycles. Below a critical noise intensity, the random area enclosed 
by hysteresis cycles is concentrated near the deterministic area, which is different for 
F*\ , small and large driving amplitude. Above this threshold, the area of typical hysteresis 

cycles depends, to leading order, only on the noise intensity. In all three regimes, we 
derive mathematically rigorous estimates for expectation, variance, and the probabil- 



ity of deviations of the hysteresis area from its typical value. 
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1 Introduction 

O 

r S^ • For a long time, hysteresis was considered as a purely static phenomenon. As a conse- 

quence, it has been modeled by various integral operators relating the "output" of the 
system to its "input" , for operators not depending on the speed of variation of the input 

> 



(see for instance [May] and [MNZ] for reviews) 

This situation changed drastically a decade ago, when Rao and coauthors published 
a numerical study of the effect of the input's frequency on shape and area of hysteresis 



cycles | RKP |. They proposed in particular that the area A of a hysteresis cycle, which 



measures the energy dissipation per period, should obey a scaling law of the form 

Ac±A a £p (1.1) 

for small amplitude A and frequency e of the periodic input (e.g. the magnetic field), 
and some model-dependent exponents a and (3. This work triggered a substantial amount 
of numerical, experimental and theoretical studies, trying to establish the validity of the 
scaling law (|l.l| ) for various systems, a problem which has become known as dynamical 
hysteresis. 



The first model investigated in [RKP] is a Langevin partial differential equation for the 
spatially extended, A^-component order parameter (e.g. the magnetization), in a ($ 2 ) 2 - 
potential with 0(A r )-symmetry, in the limit A^ — ► oo. Their numerical experiments 
suggested that ( |1.1| ) holds with a ~ 2/3 and j3 — 1/3. Various theoretical arguments 



[DT, SD, ZZ| indicate that the scaling law should be valid, but with a = j3 = 1/2. 



The second model considered in | RKP|| is an Ising model with Monte-Carlo dynamics. 
Here the situation is not so clear. Different numerical simulations (for instance |LP| , [AC] , 
ZZLJI ) suggested scaling laws with widely different exponents. More careful simulations 
| SRN|| , however, showed that the behaviour of hysteresis cycles depends in a complicated 
way on the mechanism of magnetization reversal, and no universal scaling law of the 



form (1.1) should be expected. Rigorous results on hysteresis in the Ising model are only 
available for discontinuous reversal (quenching) of the field [ JSS| ] . 

A third kind of models for which scaling laws of hysteresis cycles have been investigated 
belong to the mean field class, and include the Curie-Weiss model. A one-dimensional 
deterministic equation modeling a bistable laser, and being equivalent to the equation 
of motion of an overdamped particle in a periodically forced double-well potential, was 



considered in [JGRM]. The area of hysteresis cycles was shown to obey the scaling law 



A ^ A + e 2/3 



ri.2 



for sufficiently large driving amplitude. A similar equation governs the dynamics of the 
magnetization in the Curie-Weiss model, in the limit of infinite system size. This equation 
was examined in [fTOfl , where it was shown that the behaviour changes drastically when 
the amplitude of the forcing crosses a threshold, a phenomenon they termed "dynamic 
phase transition" . 



As pointed out in |Rao| , the difference between the scaling laws (|1.1|) and (|1.2|) can be 
attributed to the existence of a potential barrier for the one-dimensional order parameter, 
which is absent in higher dimensions. The deterministic equation, however, neglects both 
thermal fluctuations and the finite system size, whose effects may be modeled by an 
additive white noise (see for instance [ Maijl ) . Noise, however, may help to overcome the 
potential barrier and change the scaling law. 

The aim of the present work is to give a rigorous characterization of the effect of 
additive white noise on scaling properties of hysteresis cycles. For definiteness, we shall 
consider the case of a Ginzburg-Landau potential, i. e., the stochastic differential equation 



dx.c 



d_ 
dx 



1 4 

4 s 



1 2 

-x. 

2 s 



X(es)x s 



ds + adW s , 



where W s is a standard Brownian motion, and 



\(es) = -Acos(2ires), A>0. 



(1.3) 



(1.4) 



However, our results depend only on certain qualitative features of the bifurcation diagram 



and the proofs carry over to a more general setup as in [BG1, BG2]. 
In the deterministic case a 

• for A < A, 

and only if |A| < A c ), solutions of (|1.3|) are attracted by hysteresis cycles (one for each 



0, it is known JFC], pGRMj , pK| that 
+ 0{e) (where A c = 2/(3y3) is such that the potential has two wells if 



potential well) enclosing an area of order e, and with nonzero mean; 
• for A > A c + 0(e), solutions are attracted by a hysteresis cycle enclosing an area of 
order Aq+e 2 < 3 (A — Ac) 1 ' 3 , where the static hysteresis area Ao is a constant, depending 
only on the geometry of the equilibrium branches. 

For positive cr, the area A enclosed by a trajectory during one period is a random 
variable, depending on the realization W s (u) of the Brownian motion. Our aim is to 
characterize the distribution of A as a function of the parameters e, a and oq = A — A c . 
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FIGURE 1. The different regimes as a function of amplitude A = \ c +a,Q and noise intensity 
a, for a given value of the frequency e. 



It turns out that the distribution is usually concentrated around a deterministic reference 
value. We determine the expectation and variance of A. Furthermore, we estimate the 
behaviour of deviations of A from its reference value. 

One of the main results is the existence of a threshold value for the noise intensity 
a, depending on A and e: Below this threshold, the area is concentrated near the corre- 
sponding deterministic value, while above the threshold, it depends, to leading order, only 
on the noise intensity and is slightly smaller than Ao- 

There are thus three parameter regimes, as shown in Figure |l[ with qualitatively 
different behaviour of the area A. 

• In Case I, the small amplitude regime, the area is close to the deterministic value of 
order e. There is a further subdivision into Case la, where the distribution of A is 
close to a Gaussian with standard deviation Oyfe smaller than e, and Case Ib, where 
the distribution is more spread out (see Theorem [2]^ and Figure ||) . 

• In Case II, the large amplitude regime, the area is concentrated near the deterministic 
value of order Aq + (e^/ao) 2 ' 3 . In Case Ila, the distribution is close to a Gaussian 
with standard deviation of order a{e y /ao) 1 ' & . In Case lib, we can only show that A 
is concentrated in an interval of width (e-^/ao) 2 ' 3 (see Theorem 2.4 and Figure ^|). 

• In Case III, the large noise regime, A is likely to be close to a reference area A of 
order Aq — er 4 ' 3 , which is smaller than the static hysteresis area. This is due to the 
noise driving x over the potential barrier before it becomes minimal or vanishes. The 
deviation — <r 4 ' 3 does not depend on e or A (see Theorem 2J5 and Figure [Bj). 



Hysteresis does not only occur in ferromagnets and lasers, but also in mechanical 
systems displaying relaxation oscillations, such as the Van der Pol oscillator. Here additive 
noise can also have the effect of enabling jumps between stable states separated by a 
potential barrier [Fr|. Simple climate models can also display hysteresis, as has been 



observed for instance for the Atlantic thermohaline circulation Rah, Mo]. In these systems, 
the effect of small scale degrees of freedom is represented by additive noise. Our results 



describe quantitatively how noise may cause the system to switch to another equilibrium 
state, at an earlier time than expected from the deterministic approximation. 

We presented our results in detail in Section g. Section |3| contains a short description 
of the deterministic dynamics, while the remaining sections present the proofs for the 
various parameter regimes. 
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2 Results 

We consider the non-autonomous SDE 

dx s = F(x s ,X(es))ds + adW s , (2.1) 

where F derives from a periodically forced double-well potential and W s is a standard 
Brownian motion on some probability space (O,^ 7 , P). For definiteness, we shall consider 
the case 



Fix, A) = x - x 6 + A = — — - 

ox 
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(2.2) 



A(es) = -^cos(27res), A>0. (2.3) 

We introduce the notation F to,x ° for the law of the process {xt}t^t , starting in xq at time 
to, and use M to ' x ° to denote expectations with respect to P* '^ . Note that the stochastic 
process {xt}t^t is an inhomogeneous Markov process. 

Before turning to the precise statements of our results, let us introduce some notations. 
We shall use 

• y V z and y A z to denote the maximum or minimum, respectively, of two real numbers 
y and z. 

• If <p(t, e) and ip{t, e) are defined for small e and for t in a given interval I, we write 
i/j(t,e) x ¥>(£,£) if there exist strictly positive constants c± such that c~(p(t,e) ^ 
ijj(t, e) ^ c+<p(t, e) for all t E I and all sufficiently small e. The constants c± are un- 
derstood to be independent of t and e (and hence also independent of small quantities 
like a and, possibly, ao, which we consider as functions of e). 

• By g(u) = 0(u) we indicate that there exist 5 > and K > such that g(u) ^ Ku 
for all u 6 [0, 5], where 5 and K of course do not depend on e or on the other small 
parameters ao and a. 

• Let I be an interval. The notation li(x) is used for the indicator function, taking 
value 1 if x € I and otherwise. 

Finally, let us point out that most estimates hold for small enough e only, and often only 
for P-almost all w£$l. We will stress these facts only where confusion might arise. 

Let us first consider the deterministic case a = 0. It is convenient to introduce the 
slow time t = ss, and rewrite ( |2.1[ ) for a = as 

e^i = F(x u X(t)). (2.4) 
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Figure 2. Equilibrium branches of F (heavy curves) and periodic solutions of the deter- 
ministic equation (light curves), for A < A c (left) and A > A c (right). For A < A c + 0(e), 
the enclosed area is of order e while for A > A c + O(e), it is of order Ao + e 2 ^ 3 (A — Ac) 1 / 3 . 



We start by discussing some properties of this equation, which will be summarized in 
Theorem 2S2 below. As e goes to zero, solutions of (|2.4D are known to approach equilibrium 
branches of F, that is, solutions of F(x, A) = (see Figure |2|). Let A c = 2/(3\/3)- 

• For |A| < A c , F has three equilibrium branches A* (A) < Aq(A) < A* (A), where 
-X±(A) are stable equilibria and Aq(A) is an unstable equilibrium of the associated 
frozen system x = F(x, A). 

• At A = — A c , the branches A+(A) and Xq(X) undergo a saddle-node bifurcation, and 
X% (-A c ) = X *(-A c ) = x c := 1/V3. 

• For A < — A c , X*(A) is the only equilibrium branch. 

• A similar bifurcation occurs at A = +A C , where X*(A C 

• For A > A c , A/j_(A) is the only equilibrium branch. 

We can thus expect a qualitative difference, in the limit e — 

where F always derives from a double- well potential, and the regime A > A c , where F has 

only one equilibrium part of the time. 

Definition 2.1. Let x^ er ' £ be a periodic solution of (2.4). We say that this solution does 
not display hysteresis if there exists a continuous function A i— > X*(X) such that 



Xr 



' = X$ (A c ) 
0, between the regime A < A c , 



v per,e 

hm xK 



X*(\(t)). 



(2.5) 



// no suc/i function exists, we say that x^ er ' £ displays hysteresis. 

If A < A c , solutions starting near a stable equilibrium branch ACi(A) or X*(A) will 
remain close to that branch, and relation fl2.5| ) holds with X*(A) = Xi(A) or A* (A), 
depending on the initial condition. If A > A c , however, it turns out that 



v per,e 

hm xK 



A*(A(t)) if A(t) > A c or if A(i) > 
A*(A(t)) otherwise. 



-A c and \'(t) < 



(2.6) 



Thus the solution displays hysteresis since the instantaneous value of A alone does not 
suffice to determine the state of the system in the adiabatic limit. This so-called hysteresis 
cycle can be characterized by its area, defined as 

»l/2 



A(e) 



zp e A'(t)dt. 



(2.7) 



-1/2 



If A < A c , we have lim e -^0"4( e ) = 0, while for A > A c , 

"A. 
'-A c 



]imA(s)=A := [ XC (XX(X)-Xl(X))dX = l (2.8) 



The situation is thus relatively simple in the limit e — ► 0. Since in practice, however, the 
variation of A will not be infinitely slow, it is important to understand what happens for 
small but positive values of e. We summarize the necessary facts in the following theorem. 

Theorem 2.2 (Deterministic Case). There exist constants 71 > 70 > such that the 
following behaviour holds for sufficiently small e. 



If °o = A — A c ^ 7q£, Equation (|2.4|) has exactly two stable periodic solutions x 



pci 

t 



and x\ ev ~, and one unstable periodic solution x^ cr ' . These solutions track, respec- 
tively, the equilibrium branches X±(\(t)) and Xg(A(i)) at a distance not larger than 
O(c:|ao| _1 ' 2 A \^e), and enclose an area 

A{e) x eA. (2.9) 

All solutions which do not start on xf er ' are attracted either by x]? er ' or by x^ ei '~ . 
If ao = A — A c ^ 71 £, Equation ( |2.4D admits exactly one periodic solution x^ er . This 
solution is stable, satisfies ( |2.6D in the adiabatic limit, and encloses an area A(e) 
satisfying 

A{e)-A ^e 2 / 3 al /3 . (2.10) 



In the case where ao is of order 1, the scaling law (2.10) was first obtained in |JGRJ\1| ], 



We outline the proof of Theorem \L^ in Section y. Note that in the transition zone 
7o£ < ao < 7i£, the situation is more complicated, since more than two stable periodic 



orbits can coexist [TO, BK]. 



Let us now return to the stochastic differential equation (2.1). In slow time t = es, it 
can be written as 

dx t = -F(x t ,X(t))dt + ^ r dW t . (2.11) 

e ye 

Let us fix, say, to = — 1/2 as initial time, and some xo > as initial condition, such 
that the solution xf et of the deterministic equation ( |2.4[ ) with the same initial condition 
is attracted by x^ er ' or x^ er , respectively. We denote by xt(to) the solution of the SDE 



(2.11) with initial condition xt = xq for a given realization uj of the Brownian motion, 



and associate with it the area 

-1/2 

-1/2 



/■1/2 
A{e,<j;u) = - I x t (u>)X\t)dt. (2.12) 

7-1/2 



Note that A{e,o;uj) also depends on ao = A — A c . We do not stress this dependence 
here but consider ao as a (possibly constant) function of e. Of course, since xt(u>) is not 
periodic in general, the integral fl2.12|) does not represent the area of enclosed by a closed 
curve. However it is still physically meaningful since it describes the energy dissipation if 
x and A are thermodynamically conjugate variables. One can check that for \xq — x 



per, 



sufficiently small (but still of order one), \xf et — x^ er ' | decreases exponentially fast in 
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FIGURE 3. A sample path of Equation |2j for e = 0.001, a = 0.05 and a = -0.1, 
corresponding to the small amplitude regime. 



(t — to)/e and thus A(e, 0) still behaves like (|2.9|). The same is true for x^ er and the 



validity of ( pTT0| ) . 

Our main purpose is to characterize the distribution of the random variable A(e, a) as 
a function of the parameters e, clq = A — A c and a. The following three theorems describe 
the situation in three different parameter regimes. 

Theorem 2.3 (Case I — Small amplitude regime). Assume that ao = A — A c ^ 7o£ 

and that a ^ (|ao| Ve) 3 ' 4 . T/ien t/iere exist positive constants k, hi, h 2 , cq and C such 
that the following properties hold for sufficiently small e. 

• The probability that a sample path starting near one potential well crosses the potential 
barrier during one period is smaller than 



C K{ \ ao \ V£ f/2 /(T 2 



(2.13) 



Case la: Assume that either oq ^ — e or a ^ ^/e/|log|ao||. Then the deviation of the 
area A(e, o~) from its deterministic value A(e, 0) satisfies 



¥{\A(e,a)-A(e,0)\ ^ H} sC < 



- e"^ 2 /^ forO^H^ H x (e, ao), 



C 



-kH a /<t 2 



(2.14) 



for H > h 2 , 



where H\(e,ao) = h\y/e(\ao\ Ve) 3 ' 4 A(e/|log(|ao| Ve)|). Furthermore, under the slightly 
stronger assumption cr|loge| ^ co(|ao| Ve) 3 ' 4 , 



\K{A(e, a) - A(e, 0)}| < 0(<7 2 |log(|a o | V e)\) 
Var {A{e, a) - A(e, 0) } x a 2 e. 



(2.15) 
(2.16) 



Case lb: Assume now that oq ^ —e and a ^ -y/^/ 1 1°S I a o 1 1 ■ Then ( 2.14] ) still holds, 
and in addition, we have 



¥{\A(e,a)-A(e,0)\ ^ H} ^ -e 



- K H/(<7 2 |log|a | 



(2.17) 



for e/|log|ao| | ^ H ^ H2(e,ao) = /ii|ao| 3 ' 2 |log|ao||. Moreover, under the slightly 
stronger assumption <r|log£| ^ co(|ao| V e) 3 ' 4 , 

\E{A(e,ff) -A(e,0)}\ ^ C(cr 2 |log|a | |) (2.18) 

Var{A(e,a) - A(e,0)} ^ O(cr 4 |log|a || 2 ). (2.19) 

The proof follows as a particular case of more general results presented in Section |j. 

In Case la, the distribution of A(e, a) is close to a Gaussian centred at A(e, 0). Both 
the expectation of A(e, a) — A(e, 0) and its standard deviation are smaller than the deter- 
ministic value A{e, 0) x e. Thus one will still observe, with a high probability, an area of 
the same order as the deterministic one. 

In Case lb, the distribution of A(e, a) becomes more spread out, with a standard 
deviation possibly exceeding the deterministic value A(e, 0). Thus, although typical values 
of A(e, cr) will still be small, the probability of negative values is no longer negligible, and 
the deterministic scaling law A(e, 0) x e can no longer be observed. 

The quartic decay of the probability of deviations of order larger than 1 from the 
deterministic area is a consequence of the cubic growth of the drift term F for large |x|. In 
fact, this property holds in all parameter regimes, since it does not depend on the details 
of the dynamics near the origin. For the sake of brevity, we will not repeat this estimate 
in the other regimes. 

Note that there is a gap between H ^ Hi or H 2 and H ^ h 2 where we do not describe 
the deviations. In fact, the distribution of A{e, a) will not be unimodal. Sample paths are 
unlikely to jump from one potential well to the other one, but if they do so, then most 
likely near the instant of minimal barrier height, producing a small peak in the distribution 
for areas A(e, a) of order 1. 

Theorem 2.4 (Case II Large amplitude regime). Assume that oq = A — A c ^ jie 

and that a ^ {e^/aH) 1 ' 2 . Then there exist positive constants k, h\, h 2 , cq, Lq, Li, L 2 and 
C such that the following properties hold for sufficiently small e. 

• Let A denote the (random) value of A when xt changes sign for the first time. Then 

P{|A°| < A c - L} ^ - e --(i 3/2 v^)/- 2 (2.20) 

for -Li(e^) 2 / 3 ^ L ^ L /|log(e^)|, and 

P{|A°| > \ C + L) < 31 (0 ,ao](£) expj-4 f I (2.21) 

I cr 2 (e^) 2 /^|log(e^)| J 

forL>L 2 {e^) 2 /\ _ 

• Case Ila: Assume that a ^ (e^/ao) 5 ^ 6 . Then 

F{\A(e,a)-A(e,0)\^H}^-exp{- K —-^—-—) V H < h l£y /^. (2.22) 
Furthermore, under the slightly stronger assumption cr|loge| ^ Co(ey / ao) 5 ' 6 , 

\H.{A{e, a) - A(e, 0)} | ^ ( £^s ) ^ 

Yav{A(e,a) - A(e,0)} x a 2 ^^) 1 / 3 . (2.24) 
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FIGURE 4. A sample path of Equation |2dJ for e = 0.005, a = 0.04 and a = 0.04, 
corresponding to the large amplitude regime. A solution for a = is shown for comparison. 



• Case lib: Assume now that (ey/ao) 5 ^ 6 ^ a ^ {e^/ao) 1 / 2 . Then 

F{A(e, a) - A(e, 0) ^ -H} ^ - e -^ 3/2 /^ (2.25) 

F{A(e,a)-A(e,0) > +H) ^ -exp{-K ^£j°£ig- | ( 2 .26) 

el cr 2 1 log (e^) | J 

/or /ii(e A /ao) 2//3 |log(e v / ao)| ^ i7 ^ /i2(£ v / ^o) 1//3 |l o s( e2//3a o )l- ^ s a consequence, if 
a ^ co(£y^)/y^log£|, i/ien expectation and standard deviation of A(e,a) — A(e,0) 
are both at most of order (ey / ao) 2 ' 3 |log(ey / a u ")|. 



The proof is given in Section 

The estimates ( p. 20 ) and ( 2.21| ) show that the value of the parameter (e. g. the magnetic 
field) for which xt changes sign is most likely =t(A c + 0((e^/a~Q) 2 ' 3 )), which corresponds 
to the deterministic value. 

In Case Ha, the distribution of A(e, a) is again close to a Gaussian in a neighbourhood 
of A(e, 0). Both the expectation of A(e, a) — A{e, 0) and its standard deviation are smaller 
than the deterministic value of A(e, 0) — Aq x (e^/oq ) 2 ' 3 . Thus one will still observe, with 
a high probability, an area of the same order as the deterministic one. 

In Case lib, we can only show that A(e, cr) is likely to belong to an interval of size 
(ey / a u ") 2 ' 3 |log(ey / ao)| centred at the deterministic value, so that A(e, a) — Ao is not nec- 
essarily positive with probability close to 1. There is a gap between the estimates ( [2.25 ) 



and (2.26) outside this interval, and the trivial bound 1 inside the interval. This is due to 
the large spreading of paths during the jump. However, this result may conceivably fall 
short of being optimal. 

Theorem 2.5 (Case III — Large noise regime). Assume that either oq ^ e and a > 



(|ao| V e) 3//4 , or ao ^ e and a > (e 
area A, satisfying 



ao) 1 ' 2 . Then there exists a (deterministic) reference 



A- An 



-a 



4/3 



(2.27) 



and positive constants k, hi, /12, cq, c\, C2 and C , such that the following properties hold 
for sufficiently small e. 
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Figure 5. A sample path of Equation |2.1| for e = 0.005, a — 0.16 and a = —0.01, 
corresponding to the large noise regime. A solution for a — is shown for comparison. 

Case Ilia: Either a§ ^ e or a > a Q . Then the deviation of the area A(e, a) from 
the reference value A satisfies 

F{A(e,a 






¥{A{e, a)-A>+H}^^ e -^/^\^\) +h [QMa4/S) 



(2.28) 

( H } e -KH/(s\\0g<T\) ( 229 ) 



for ^ H ^ /iia" 2 ' 3 |logo"|. Moreover, if the noise intensity satisfies c\£ < o" 4 ' 3 /|logo"| 2 
and a 2 ' 3 1 logo" | ^ C2/|loge|, then 

E{A{e,a)-A} £ [-Ca 4/3 |loge| 2/3 , C(e V a 2 |loge|)|logo-|] (2.30) 

Var{A{e,a)-A} < C(o- 4/3 |loge| 2 / 3 ) 2 . (2.31) 



• Case Illb: oq ^ e and a ^ a . Let £q = |log(c 4 ' 3 /\/ao)|. Then 



»{A(£,a) -A < -H} ^ C e ^ H 3/2 /a 2 + 3 e _ KCT 2 /(Ev ^| logCT 

C — ■ 



I) 



(2.32) 



-,.,///(a 2 £ )_^ 1r /^ e - K a 2 / 3 H/( ev ^|log CT |) ^33) 

/or < H < /iio -2 / 3 ^- -fri addition, if a ^ ci(e v /ao) 1/ ' 2 |loge| ; i/ien 



»{.A(e,<7) - i > +#} < - e -»"/^«o; +-l [0 , h2ao) ( 



E{^(e,o-) -i} e [-Co- 4/3 |loge| 2/3 ,C(o- 2 4|loge| VeVooliog^l/^ 73 )] ( 2 - 34 ) 
Ym{A{e,a)-A] ^ C(o- 4/3 |loge| 2/3 ) 2 . (2.35) 

The value A of A ai the first time xt reaches behaves in a similar way as the area, when 
compared to a reference value A equal to A c — 0(o" 4 / 3 ). 

The proof is given in Section & 

The main feature in this parameter regime is that the noise intensity is sufficiently 
large to drive xt over the potential barrier before it reaches its minimal height or even 
vanishes. The barrier is typically crossed when |A| equals A c — 0(a A ' 3 ). 

The distribution of A(e, a) decays faster to the right of A than to the left. The 
probability that A(e,a) exceeds .4.0 is very small (unless a approaches its threshold value), 
so that it is indeed likely to observe an area that is smaller than Ao, by an amount of 
order o -4 ' 3 . 



10 



3 Deterministic case 

In this section we discuss the deterministic equation 

e-^- = x - x 3 + X(t) (3.1) 

at 

with initial condition x_i/ 2 = ^o- Recall that we are interested in the case X(t) = 
—Acos(2irt) with A = A c + ao- Since this equation has already been studied in |JGRM| , 
TO , |BK] , |BG2| , we only outline the main properties without proofs. 



3.1 The case a ^ 7 e 

The simplest situation occurs when ao is negative and of order 1. Then the three curves 
x^_(t) = X±(\(t)) and Xq(£) = Xq(A(£)) are uniformly hyperbolic equilibrium curves of the 
associated family of frozen systems x = x — x 3 + A. Thus Tihonov's theorem fTi| , |Gr[ shows 
the existence of particular solutions x t and x® tracking, respectively, x± (t) and Xq (t) at a 
distance of order e. These solutions are not necessarily periodic, but the curves x t attract 
a neighbourhood of order 1 exponentially fast. Thus the Poincare map P : x^\/2 *— ► £1/2 
maps neighbourhoods of order 1 of x±(— 1/2), respectively, to two exponentially small 
intervals containing x, ,„. This implies the existence of a unique fixed point (corresponding 
to a periodic orbit) in each interval. A similar statement is true for P _1 in a neighbourhood 
of Xq(1/2). The fact that x t is monotonous between the equilibrium branches excludes the 



existence of other periodic orbits (see, e.g. |Ber| , Proposition 4.8]). 



If ao is a small parameter, Tihonov's theorem can also be applied outside a given 
interval [—T,T], T a constant of order 1, while the dynamics in [— T, T] has to be analysed 
separately. For |t| ^ co(|ao| Ve) 1 ' 2 , cq a sufficiently large constant, one can consider the 
deviation yt = xt — x+(t), which obeys the equation 

ef t =al(t)y + bX(y,t)-e^±, (3.2) 

where 

a;(t) = l-3«(t)) 2 
b\(y,t) = -y 2 [3x\(t)+y]. 

This equation is used to show (see [ p3G2j , Section 4.1]) that 



(3.3) 



yt x ^r for -T ^ t < -co(|ao| V e) 1 ' 2 . (3.4) 



In the case ao ^ — £/7o for some small enough 70 > 0, one also obtains from ( |3.2| ) that for 
1*1 < cov / |ao|, 

yt = -tC 1 (t) + C 2 (t) with d(i)x^, C 2 (i)x-^-7o, (3.5) 

«ol l«o| ' 

which implies that yt becomes negative at a time of order e|ao| -1 ' 2 . If — e/70 ^ ao ^ 7o£) 
the dynamics for \t\ ^ coy^ is analysed by the change of variables 

t = c\fes, x = —={\-\ z J, (3.6) 

v 3 ^ c / 
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where 2tt 2 Ac a = 1. Then z obeys a perturbation of order \fe of the Riccati equation 

— = s 2 — z 2 — 5, where 5 = V3c 2 — < y/3 c 2 70, (3.7) 

as e 

which can be used to show that for 70 small enough, z s x 1 for s of order 1. It follows 
that for all ao ^ 7o£, 

xt - -j= x (|ao| V £) 1 ' 2 for |t| < co(|o | V £) 1 ' 2 . (3.8) 

Finally, one obtains as before that 

Vt - - ||r for cd(|ao| V £) 1 ' 2 < * < T. (3.9) 



Hence there is a solution of ( |3.lD tracking x* (i) at a distance of order er/(|t| V y|oo| V y/k~) 
(if ao > 0, x+(t) does not exist during a small time interval, but this gap is too small 
for xt to slip through) . Later we will use the fact that the linearization of F around this 
solution satisfies 

q _ 

a(t):= — [x-x 3 -Acos(27rt)] x -(\t\ V V|ao| V Vi). (3.10) 



dx 

Furthermore, we will need that fact that \a'(t)\ is bounded above by a constant independent 
of £ and ao . This can be shown by using the relation 

a'(t) = -6x tT -x t = -6x t -F(x t , \(t)). (3.11) 

at £ 

The cases \t\ ^ (|oq| V e) 1 ' 2 and ao ^ — e/70 can be treated by expanding F around the 



equilibrium branch x* + (t) that xt is tracking, and using the estimates (3.4), (3.5), ( |3.9D for 



xt — x\(t). The remaining case can be treated by considering (|3.7| ) directly. 
In addition, these estimates show that 



i.1/2 

A(e, 0) = 2ttA / x t sin(2vrt) at x eA. (3.12) 

7-1/2 



Similar properties hold for solutions tracking x* (i) and Xq{€), and the above arguments 
on the Poincare map can be repeated to show the existence of two stable and one unstable 
periodic orbit. 

3.2 The case a ^ 716 

Let t c be the solution of ^4cos(27r£ c ) = A c in [0, 1/4]. The equilibrium branches x* + and Xq 
bifurcate at the point (— t c , x c ), where x c = l/\/3. The translation 



yields the equation 



-t c + s, x = x c + y (3.13) 



e^ = ls(s)-V3y 2 -y 3 , (3.14) 
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where 



fi(s) 



A c - 
A c (l 



^lcos(27r(-t c + s)) 

- cos(27rs)) - ^J A 2 - A 2 sin(27rs) 



2irJ2a \ c + a^s + 2-k z X c s z + 0(s 6 



x: 



(3.15) 



"6c "T" S) Xq at a 



As before, one shows that y s tracks the equilibrium branch y+(s 

distance scaling like e/\s\ for s ^ —e 2 ' 3 a . For larger times, we use the scaling 



2/3-1/6 



ce 'a 



u, 



y/3t 



-s^a^z 



(3.16) 



which yields, for an appropriate choice of c, a perturbation of order e 1 ' 3 a of the Riccati 
equation 



dz 
du 



fi(u) - z 2 , 



where ji{u) = -u + O(e 2/3 a 2/3 u 2 ). 



(3.17) 



For sufficiently large 71 (recall that ao ^ Tie); one can show that z u reaches a value of order 
— 1 in a time of order 1, and thus y s reaches order — e ' 3 a for some s of order e 2 ' 3 a . 
Finally, the fact that the right-hand side of ( 3.1 4| ) is smaller than — y 2 for sufficiently small 



s and y can be used to show that y s reaches values of order —1 after another time of order 



r2/3„-V6 



For larger times, y s is quickly attracted by the lower stable equilibrium branch. 



Later we will use the fact that the linearization around y s satisfies 







a(-t c + s) := — [/x(s) - \/3y 



dy 



y=ys 



(|s| VaJ /4 ^Ve 1/3 aJ /6 ) (3.18) 



-1/6 



for s ^ Co£ 2 / 3 a, 

constant times (ao/e) 1 / 3 for t < — t c + O(e 2//3 ao 



Furthermore, one can check that |a'(i)| is bounded above by a 

. It is easy to see that the solution we 



1/6 



constructed encloses an area satisfying A — Aq >c e 



2/3^/3. 



where the main contribution 



comes from the delayed jump from a neighbourhood of x\(t) to a neighbourhood of #*(£). 
Since a;+(t) is the only equilibrium branch near to = —1/2, the Poincare map contracts 
any interval of order 1 containing x\{t) to an exponentially small neighbourhood of £1/2, 
which implies the existence of a unique periodic orbit. 



4 The random motion near stable equilibrium branches 

We consider now the SDE 

dx t = -\x t - x 3 - Acos(2irt)] dt + -^= dW t (4.1) 



with a given (deterministic) initial condition xt Q = xq . Let xf et denote the solution of the 



deterministic equation (|3.l| ) with the same initial condition. We will start by investigating 



the difference Xt 
difference. 



x 



dot 



, and then derive some properties of the area delimited by this 
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4.1 Noise-induced deviations from the deterministic solution 

The difference yt = xt — xf ct obeys the SDE 

dyt = -[a(t)y t + b{y t ,t)]dt + ^ F dW t , y ta = 0, (4.2) 

£ ye 

where 

„det\2 



a(t) = 1 - 3(xt 
b(y,t) = -y 2 [3xt et + y]. 



(4.3) 



In this section, we are interested in situations where x^ et is attracting up to time t, that 
is, we assume that a(s) < for to ^ s ^ t. Results from the previous section (see fl3.10| ) 
and ( |3. 18[) ) show that this is true when the following condition is satisfied. 

Assumption 4.1 (Stable case). Assume 

• either oo ^ 70 £ and t arbitrary 

• or ao > 7i£, to ^ —1/2 and t ^ t* ■= — t c + coe 2 ^a^ 

For the sake of brevity, we will refer to these assumptions as stable case. 

If we were to omit the nonlinear term b(y,t) in ( |4.2| ), the solution yt of the equation 
would be normally distributed with mean zero and variance 

a 2 r l /"' 

v (t) = — e 2a ^' s ^ £ ds, where a(t, s) = a(u) du. (4.4) 

£ Jto J s 

It is straightforward to show that the function 

C ( t ) := -J— e Mt,to)/e + 1 [* e 2a(t,s)/e ds (45) 

2|o(*o)| e Jto 



satisfies, in both cases summarized in Assumption 4.1 



m x -i-. (4.6) 

2\a{t)\ 



Thus in the linear case, the standard deviation of yt is smaller than 0"a/£(£). The following 
result applies to the whole path {y s }t ^s^t of the nonlinear equation fl4.2j) . 

Proposition 4.2. iet 

C(t) = sup COO- (4.7) 

In the stable case, there exists a constant ho (ho x 1) such that for all h ^ /ioC(t) -3 , 

'i,s<^ > A M^ + 2 ) -B50 - »w - »(«)}■ 

(4.8) 



Proof: The case ao ^ is the one considered in [ BG2] , Theorem 2.6], and the other cases 



can be proved in exactly the same way. □ 
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If we do not care for the precise value of the exponent in ( [L8|) , an obvious modification 
in the proof yields the bound 

pft>,o J sup J^fL >h \ < cf 1 -^- + l) e- Kh2 /°\ (4.9) 

Uo<s<* VC(s) J V £ J 

for all h ^ /ioC(^) -3 > where C and k are positive constants. This estimate shows that in 
the time interval [to, t], the typical spreading of paths is of order <jy/£(s) for a <C CC0~ 2 ■ 
It allows to bound the probability of deviations up to order £(t) _1 . On the other hand, the 
special (cubic) form of the drift term in Equation ( |4.2| ) allows for a bound on deviations 
of order larger than 1. 

Proposition 4.3. There exist constants Lq,C,k > such that for all L ^ Lq and all 
Vo ^ Lq/2, 

P*o»«>( sup y s > l\ ^ c( ^— ^ + 1 J e"^ 4 /' 72 . (4.10) 

Proof: First note that we are working in slow time. The estimate is classical for t—to ^ e, 
and starting from there, ( [4.10 ) can be obtained by considering a partition of the interval 



[to,t] with spacing proportional to e/(t — t ). D 

Remark 4.4. Note that the preceding proposition holds for all ao, but Lq may depend 
on the amplitude A. 

The following proposition gives bounds on the moments of yt- These bounds hold 



whenever the estimates (|4.9|) and (|4.10|) are satisfied, and we do not need to assume that 
a{s) < holds for all s. 

Proposition 4.5. Fix t > to such that t — to * s a t most of order 1 and assume that there 
exists an ho(e,t) > such that ( |4.9[ ) holds for all h ^ ho(e,t). Then there exist constants 
K, M > such that 

(i-i __ \k\ 
i + ! °^ £| J , fceN, (4.ii) 

whenever a satisfies Kay/\log(ea 2 )\ ^ ho(e, t). 

Proof: We will only prove the case k = 1, as the general case follows along the same 
lines. Let 7 = Ka^J\\oge\ for some constant K > to be chosen later and set h = ho(e, t). 
Note that, under our condition on a, we may assume 7 < h. Let Co = mi "se[to,t] C( s )- We 
write the expectation of yf as E\ + E, where 

E l= ^4y 2 t l { _U^ <7 )} and E = E*°'°L 2 lr j^L^lj- (4.12) 

The first term E\ can be estimated trivially, namely by E\ ^ 7 2 £(t). To estimate the 
second term E, we employ integration by parts, thereby obtaining 

E^((t)[ 2zF* (1 '°{ sup pj_ > 7 V z 1 dg. (4.13) 

JO Uo<s<t VC(s) J 
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We now split the integral at 7, h and Lq/^/Qq and estimate the resulting terms separately. 
By (ED, 



£2 = m I 

Jo 



Lo/V& 



2z¥ to '° 



<c(t)C ['LJ° + , 



sup — s ^ 7 V z \ dz 

2 , A - K7 2 /" 2 . ±0 p -^ 2 /^ 2 



7 H I e 



i0 - 



Estimating the remaining part of the integral with the help of ( 4,1C| ) 



E 3 = ((t)l 2zF t °'°\ sup -^=> 7 v4<lza(i)C| — + l|r 7 i 7 <' 



t-to 



cr 



(4.14) 



-kZ,4/<t 2 



(4.15) 



follows. Since the expectation of y\ is bounded above by E\ + E2 + £"3, ( 4.11 ) follows from 
the fact that we can choose K large enough to bound all three terms by some constant 
times <7 2 C(£)|loge|. □ 

In the stable case, which is our major concern in this section, the previous bound can 
be improved as follows. 

Corollary 4.6. Fix t such that t — to is at most of order 1. In the stable case, there exist 
constants c\ > and M > such that, i/(r|loge|C(t) 3 ' 2 ^ c\, then 



E to '°{\y t \ 2k } ^k\M k a 2k ((t) k , fceN. 



(4.16) 



Proof: Let us again focus on k = 1. Estimate ( |4. 16 ) is obtained in the same way 
as ( 4.11 ), the only difference lying in a more elaborate bound on E\. We use the integral 
representation 



/// "= - / <' 
e -/to 



a{t,s)/e 



b(y s ,s)ds + 



a 



3 «(M)/e dW 



(4.17) 



to 



of the SDE (\L$) (for y to = 0), thereby obtaining 



E 1 ^ 2E 



*o,o 



^a{t,s)/e 



b(y s ,s)ds) 1 



/o 



'! 



suPi «: s! S 



/CO) 



=57 } 



+ 2- 



r 2 r* 



to 



3 MM)/e ds . 



(4.18) 



The second term on the right-hand side is bounded above by 2o~ 2 ( i (i). The first one can 
be estimated by bounding b(y s , s) uniformly in s, with the help of the estimate \b(y, s)\ ^ 
Mo{y 2 + |y| 3 ), valid for all s, c.f. (|4.3| ). The remaining integral behaves like eC{t). Thus 
we obtain 



Eh. ^ constC(t) 2 ( 1 2 C(t)+ 1 3 C(t) 3/2 ) 2 + 2a 2 at), 



(4.19) 



while E can be estimated as before. Again choosing 7 = Ka^J\k>g£\ for K large, yields 
Estimate ( pUP - D 
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4.2 Noise-induced deviations from the deterministic area 



Let us now examine the behaviour of the surface delimited by the process yt- We want to 
control the process 



Y t 



y s X'(s)ds, 



(4.20) 



to 



which measures the deviation of the area A(e, a) enclosed by xt from the one enclosed by 
. Using the representation ( |4.17| ) of yt, we obtain, by a version of Fubini's theorem, 



det 



that 



where 



Y t 



g(t,s)b(y s ,s)ds + ^ g(t,s)dW i 



to 



V~e 



u, 



g(t,s) 



a a(u,a)/e \i 



A (u) du. 



In particular, the term 



i? = 4= / g(t,s)dw : 

V £ J to 



to 

is a Gaussian random variable with mean zero and variance 

r 2 ft 



(4.21) 
(4.22) 
(4.23) 
(4.24) 

(4.25) 
(4.26) 
(4.27) 



Proposition 4.7. Under the assumptions of Proposition \^.3j , there is a constant M\ > 
such that for all k ^ 1, 



Var(Y 4 °) = - [ g(t,s) 2 ds. 
£ J to 



In the sequel, we will use the following abbreviations: 

i r l 

Ti(t,t ) = - / \g(t,s)\((s) l ds, *€{1,§}, 

£ J to 

1 /■' 

T(Mo) = - / \g{t,s)\ 2 ds, 

£ Jto 
A(t,t )= [ \\'(s)\ds. 

■Jto 
In addition, we denote by (2k — 1)!! the product ]X =1 (2* — !)• 



E t0 >°{(Y t °) 2k } = (2k - 1)!! (a 2 eV(t,t )) k 
E to,0|| yt 0|2fc-i| = k i2 k - 1 Jl(a 2 eT(t,t )f k - 1]/2 

K to '°{\Y t -Y t °\ k } ^k\M^(a 2 T 1 (t,t )) k 



(4.28) 
(4.29) 



1 + 



I log el 



l + k\a 2 



2kf T 3/2(t,to) 

Ti(t,t ) 



21: 



k\ 



1 + 



Ijogej 
k\ 



1/2 



(4.30) 
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Proof: Since ( 4.28 ) and ( |4.29 ) are an immediate consequence of the fact that Y® is 
Gaussian with variance ( f4.24| ), we only need to prove (4.30). We restrict our attention to 
the case k = 2 as the case k even follows by an obvious adaptation and the case k odd is 
obtained from the case k even by an application of Schwarz' inequality. First note that 



M 2 /"* /** 

E'°> {(y t -y t ) 2 }^ / / \g(t, U )\\g(t,v)\E t 0' {(yl + \y u f)(y 2 v + \y v f)} dvdu. 



(4.31) 



Estimating the expectation of the product by Holder's inequality and Proposition |4.5| . 
(|430|) follows. D 



Remark 4.8. In the stable case, under the assumptions of Corollary 4.6, the bound ( f4,30D 
simplifies to 



rQ\k~ 



2fcA/.Nfc\l/2 



E i0 '°{ \Y t -Y t °f}^k\ Ml{a 2 Y x {t, *„)) (l + k\ a 2k ((t) k ) 



(4.32) 



The following proposition gives bounds on the probability that the deviation of the 
area from the corresponding area in the deterministic case is large. 

Proposition 4.9. 

• Assume that there exists an ho(s,t) > such that ( |4.9| ) holds for all h ^ ho(e,t). 
Then there exist constants hi, k, C > such that for any p G (0, 1), 



p to,0||y t | > H j ^ exp 



2 tt2 



(1-P) 2 H 

2T(t,t )a 2 e 



+ C 



t-t 



+ 1 exp 



K pH 

'^ri(t,t ) 



(4.33) 



whenever pH < hiFi(t, t )(h (e,t) 2 A((ty l ). 
In addition, 

¥ to,0^ Yt \ >H }^ c (tll°+ , 



whenever H ^ LoA(t, to). 



h^&l <-) 



Proof: Consider first the case pH small. By (4.21), we have for any p £ (0, 1) 

F t0 ' {\Y t \>H}^F t0 > {\Y t °\ >(l-p)H}+F t0 '°\- [ \g(t, s)\\b(y s , s)\ ds > p h\. 

18 Jtn J 



(4.35) 



The first term on the right-hand side immediately yields the first term in ( 4.33[ ) due to 
the Gaussian nature of Y t °. We denote pH/M by Q. For q € (0, 1), the second term can 
be bounded by 



P*o,o^ sup Jy^ ., 



to 



^i a/C(s) \Fl(t,to) 



1/2 



0*0)0 



\Vs\ 



-HF*"'^ sup ~y==> . .. 

to^s^t VC(S) \ i 3/2 



i-^ y /3 i 

3/2 (Mo) y J 



(4.36) 
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We choose q by 



Ti(Mo) 



ho(e,t)/\C(t)-y 2 ^3/2{t,to) 



(4.37) 



and estimate both summands in ( 4.36| ) by (|4.9|) . Note that the first summand dominates 
the second one by our choice of g, since we assumed pH ^ hiTi(t,to)(ho(e,t) 2 A C(i) _1 ). 



Thus we obtain the bound ( 4.33 ). 

For H large, we employ the trivial bound \Y t \ ^ (sup t ^ s ^ t |y s |)A(£, to) together with 



□ 



(4.38) 



Estimate fl4.10|) , thereby obtaining (|4-34| ) . 
Choosing p = p(H, e) in ( [4.33 ) by 

v _ h ri(Mo) 

(l-p) 2 2«er(t,to) 
yields the following corollary. 
Corollary 4.10. There exist constants h\,K,C > such that 

P ,,o {|1 . ll>H} , (^£ + 1 ) exp {_yi!^} (4 , 9) 

for all H satisfying (1 — p) 2 H 2 ^ 2KhieT(t,to)(ho(e,t) 2 A C(t) -1 ). i^ere p zs defined by 



( 4.38 ). Furthermore, whenever H ^ const eT(t, to) /T\(t, to), then 1 — p is bounded away 
from zero. 

In order to complete the proof of Theorem [LQ, we have to control the function g(t, s), 
defined in ( [4.22 ). This task is simplified by using the following lemma. 



Lemma 4.11. Assume that \a'(u)\ ^ a\{e) anda(u) ^ —cy/eai(e) for allu in an interval 
[s,t], where c > is independent of ao and e. Then there exists a constant d > 0, 
independent of s, t, c, such that 



g(t,s) 



\a(s)\ 



X'(s) + O 



a(s 



(4.40) 



whenever t — s ^ de/\a(s)\. 



Proof of Theorem ^]3|. First note that Proposition [O] establishes the bound ( |2.13|) 
on the probability of a sample path crossing the potential barrier. From ( p. 10 ), ( [4.6| ) and 
the preceding lemma, one easily obtains 



r 1 (i,-i) = 0(|log(|ao|V£)|) and r(|,-|)xl 



(4.41) 



in the stable case, while A(i, — o) x 1 is trivial. Now the preceding results imply the 



stated estimates. 



□ 



19 



5 The large noise regime 

In this section, we consider those parameter regimes in which the noise intensity a is large 
enough to allow for transitions from one potential well to the other one, with a probability 
close to 1. Depending on the amplitude, there are three cases to consider: 

• ao ^ — e and a ^ |ao| 3 ' 4 ; 

• |oo| ^ £ and a ^ e 3 / 4 ; 

• ao ^ £ and a ^ {e^/oq) 1 ' 2 . 

Actually, we will need to assume that a ^ K|ao| 3 ' 4 , o~ ^ Ke 3 ' A or a ^ if (E-^/ao) 1 ' 2 , 
respectively, for some large constant K, but in order not to overburden notations, we will 
assume that K = 1 is a possible choice. 

Recall that in the deterministic case, transitions are impossible if ao ^ 7o£, and occur 
only after time —t c + co£ 2 ' 3 a if ao ^ 7i£. It turns out that under the above conditions 
on a, transitions are likely to occur some time before the potential barrier reaches its 
minimal height or even vanishes. For brevity, we shall only discuss the case |ao| ^ e in 
detail, but the other cases can be investigated similarly (since transitions occur early, they 
are not influenced by the details of the bifurcation or avoided bifurcation). 

5.1 The transition time 

We assume |ao| ^ £ unless stated otherwise. By symmetry, we may restrict our attention 



to a half-period, say t € [—1/4, 1/4]. Let xt be the solution of the SDE (|2.1l|) starting at 
time to = —1/4 in the upper well, i.e., near ori(to). We define the transition time as the 
stopping time 

r° = infjs > t : x s < 0} £ (to,oo], (5.1) 

when x s crosses the t-axis for the first time. The choice of x s = is purely for convenience, 
and the qualitative behaviour of r° remains the same if is replaced by any level between 
— x c — 5 and x c + 5 as long as 5 > is chosen in such a way that f(x,t) ^ holds for all 
|x| ^ x c + 5 and all t in question. The following result characterizes the distribution of r°. 

Proposition 5.1. There exist constants C , c±, oi, K > such that 

• fort <t < -cict 2 / 3 , 

P .»,«. {T o <(Hc (^£ + 1 ) exp {__JL_} ; (5 . 2) 

• for -ci<7 2 / 3 + c 2 £ ^ t < ci<r 2 / 3 , 

F t0 ' X0 \T°>t}^-exp\--^—-f |a(s)|d4+e~ K/<T2 . (5.3) 

2 I | log a | £ J_ ClCT 2/3 J 

det,H 



Recall that a(s) is the linearization of the drift term along x s ' as defined in ( 3.10 ), and 
((t) is defined in (|4.7|) . 

Proof: The first part is a direct consequence of (|4.9|) with h = hi((t)~ 3 ' 2 , where hi is 
chosen sufficiently small that the relation \x s — x^ et \ ^ hi/Q{s) for all s G [to,t] implies 
that i, > for these s. 
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The second part is an application of Theorem 2.7 in | BG2| (with h = const cr\log e] 1 



1/2^ 



Note that the theorem naturally extends to the case < ao ^ e. In fact, the integrand in 
(|5.3| ) should be the curvature of the potential at the deterministic solution tracking the 
saddle x^(t), but the curvature behaves like |a(t)|, compare |BG2, Proposition 4.3]. □ 



The condition a ^ e 3 ' 4 implies that C(i) _1 x l a (*)l x 1*1 f° r * ^ ~~ c\a 2 ' 3 , and thus the 
exponent in (|5.2| ) scales like |i| 3 /<7 2 . The integral in ( |5.3[) behaves like a 2 ' 3 (t + c±a 2 ' 3 ). 

Proposition JO| shows that the transition is likely to occur close to time t\ = —c\a 2 ' 3 , 
which satisfies A(ii) + A c x <r 4 / 3 . Therefore, we should compare the area A(e,a) to a 
reference area A given by 

l A= f tl x ^+(-\'(s))ds+ [ t2 x^>-(-\'(s))ds, (5.4) 

1 Jt Jt\ 

where to = —1/4, t\ = —c\a 2 ' 3 , £2 = 1/4, and we denote by x s c ' the deterministic 
solution starting in xq, which tracks x+(s), and by x s e ~ a deterministic solution tracking 
x*_ (s). It is easy to check that 

A-Ao^-a^ 3 . (5.5) 

(This relation does not depend on the initial conditions of x s e ' , as long as they are 
sufficiently close to x\(to) or x*_ (ii), respectively). 



Remark 5.2. Proposition ^lj also holds for ao < — e, with the same exponents. 

In the case ao > e, C(*) _1 behaves like |a(i)|, given by ( |3.18| ). The bound flO) holds for 
to < t < t\ = -t c -c\(a 2 / 3 A<7 4 / 3 aQ 1/2 ), with the exponent replaced by -koq \t+t c \ 3 / 2 /a 2 
if a ^ a and by — n\t + £ c | 3 /o" 2 if <r ^ a . Note that in both cases, \(t\) + A c x u 4 ' 3 . 

The bound (5.3) holds for ti + C2£ ^ t ^ —t c + cq£ 2 ' 3 a , with an exponent of the 



same order as in the other cases, namely a 2 ' 3 (t — t\) / (e\log a\) . The behaviour for larger 



t will be discussed in Proposition |6,l| below 



5.2 Deviations from the reference area 

Our aim is to characterize the deviations of the random variable A(e, o~) from its determin- 
istic reference value A over one half-period. We focus again on the case |ao| ^ £• With a 
slight abuse of notation, we can write this deviation as ^ (A(e, a) — A) = Yj£ + Y t ~ , where 

Y+:= f^^-xf^i-X'^ds, Y t -:= f\x s - xf t -)(-A'( S ))d«. (5.6) 

J to J t\ 

We will estimate separately the probability that each of these terms is larger than H or 
smaller than —H. To do so, we need a preparatory result allowing to extend the estimate 
(|4.9|) to larger values of h. 



Proposition 5.3. Define the stopping time 

t = inf{t G [t , h] : x t ^ xf t,+ - /loCW" 1 } 6 [to, *l] U {oo}, (5.7) 

where the constant ho is taken from, Proposition \4-^- Then 



I _ dct,+ i \ (~i 

p*o,xoj sup \xs-xs_ | >h \^^ e -. h y*i (58) 

Jo^sCt VC(s) J £ 



for some C, k > 0, all t £ [to,*l] an d & H h > 0. 
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Proof: The fact that the drift term F has a negative second derivative with respect to x 
for all x > implies that x s is unlikely to exit the strip of width h^jQ(s) through its upper 
boundary, as was proved for negative oq in [BG2, Proposition 4.5]. We also know by (|4.9|) 
that x s is unlikely to exit the strip through its lower boundary if a <C h ^ /ioC(i)~ 3 • 
The stopping time r has been defined in such a way that x s cannot leave a strip of larger 
width before time r. □ 

Note that by decreasing ho if necessary, we can arrange for r < tq. We are now able 
to estimate deviations of Y^. 

Proposition 5.4. There exist constants C , k, h\ > such that 



C 



-nH 3 / 2 /a 2 



C 



-Ktf/0 2 |log<x[) 



P *o,s jy+ < -H) < -e 

F t ,x ( Y + > +if | ^ e - K fl 2 /(" 2 ^) + H, 

/orO ^ H ^ /iicr 2 / 3 |logcr|. 
Proof: We decompose 

+ E i ^'|l {T<il} P^|£ 1 (x s -xf t ' + )(-A / ( S ))d s < 



(5.9) 
(5.10) 



(5.11) 



^ 



where r is defined in fl5,7|) . The first term on the right-hand side can be estimated as 



in Proposition 4.9, as there is no need to distinguish positive and negative deviations for 



this term. However, Proposition |5.3| allows us to obtain bounds valid on a larger domain 
of H. Note that Proposition 4.6 remains valid when Yt is replaced by Y rA t- This is a 



consequence of the r,(t,to) being monotone functions of t € [£o>£i] an d a slightly more 
elaborate estimate showing that sup 4o<s<f |y s | obeys the same bound as was used for |Y 4 °| 
in ( 4.35| ). Thus we obtain the estimate 



P i0 ^{|Y+ tl |>tf}^exp 



(1-p) 2 H 2 \ C 
2T(t 1 ,t )a^ej + e 6XP 






(5.12) 



valid for pi? ^ const Ti(ti,to)/((ti). An application of Lemma [4. 1 1| shows that r(ix,£o) x 
1, ri(ti,to) x | log ex | , and we already know that C(£i) x 0"~ 2 • Choosing pxl provides 



an estimate of the form (5.10). The second term on the right-hand side of (5.11) can be 
estimated, using the monotonicity of A in [io^i]; by the relation 



[j\ Xs - x ^+)(-X\ s ))ds<-H} 

^ p t ' Xt | 



sup \x s 



>L\ + 



}+P r ^{A(r)-A(£i)>tf/L}. 



(5.13) 



The first term on the right-hand side can be estimated by Proposition 4.3, provided L is 
larger than some constant of order 1. It decreases like e _cons *' CT je. Using the fact that 
A(r) — \{t\) x t 2 — t\ for r not too close to to (note that the contribution of r close to to 



is even smaller) and (|5.2| ) of Proposition 5.1] , we obtain that 

C 



E*o-*o/l {T<tl} P T ^-{A(r) - A(ti) > H/L}} < <^ e ~ const h*/*/** _ (5>M) 
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This last term is easily seen to dominate all others, so that (5J)) is proved. 

To estimate deviations in the positive direction, we split terms as in ( 5.11| ). The first 
term can also be bounded by ( |5.12|) . Using the fact (compare |BG2 , Proposition 4.5]) that 

it only remains to estimate 



l' dct.+ \ j~, 

■ sup Xa ~^ >^Hf2\^-^ H/(7 \ (5-15) 



r pt\ _ aet,+ a 

™\ (* s -^ et ' + )(-A'( S ))ds>#/2, sup Xs ~^ ^VW/2\ 



«i 



^P T ' XT {|r| 3/2 > a + constVH}. (5.16) 

The expectation of this term also decreases like q- kH / u js by Proposition 5.1. Taking 



( 5.12 ) and ( 5.15J ) into account, we have proved ( |5.10| ). □ 



The term Y t2 can be controlled in a similar way: 

Proposition 5.5. There exist constants C, k, hi, h% > such that for all xt ± £ [—L,L] 
(L x 1), and all H ^ hs, 

pil,x tl (y- < _tf\ <; £ e -*H/a 2 + 3 e - K(J 4 / 3 /(e|log<7|) + e ~ K H 2 /(a 2 e) (g^ 

Proof: The proof being similar to the one of the previous proposition, we only outline 
the main steps. Introduce a stopping time r c = inf{s G [£1,^2] : x s ^ — x c — 5} for some 
small 5 > 0, cf. the comment on the definition of r° in the beginning of the subsection. 
We first need to control the behaviour of 

/■t c A0 rr c At 2 

Y- At2 = / (x s -4 et <-)(-A'( S ))d S + / (* s -4 et '-)(-A'( S ))d5. (5.19) 

Jti Jt c A0 

Observe that since x s > x s ' for s ^ r c , the first term on the right-hand side is positive, 
while the second one is negative or zero. First note that if x s is bounded above by L x 1, 
then Y~ At cannot exceed a value of order a A < 3 . Deviations of Y~ At2 in the positive 
direction can be bounded using a decomposition similar to ( pM3 ) and applying (5.3) for 
r c instead of r. We find 

^^{Y- M2 > H) ^e"^ 2 +h [0Ma4/3) (H)e^ H ^ lo ^\ (5.20) 

valid for H > 0(ea 2 / 3 ). 

Deviations of Y~ At2 in the negative direction can only be caused by the second term on 
the right-hand side of (5.19J). However, there is no small lower bound for that term. The 



reason is that transitions to x s e ' are only probable in the window t 6 [— ciO" 2 ' 3 ,cict 2 ' 3 ]. 
If this opportunity is missed, which happens with a probability of order e~ Kcr /y £ ^°s cr \)^ 

TcAt 2 



then x s keeps tracking x s e ' i ~ and Y T M may reach negative values of order 1 
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To complete the proof, we need to show that on {r c < £2} 



1-2 



\ Xs - x f^-)(-\\ s ))ds 



>tfUe-^ 2 /( CT2£ )+-e-'^ 2 



(5.21) 



det,r c 



Let x s ' c be the deterministic solution starting in x Tc at time r c . This solution is attrac- 



D 



tive, and thus ( 5.2 1] ) holds with x s e ,Tc instead of x s e ~ as a consequence of Proposition 4.S . 
But the distance between x s e ' Tc and x s e ~ decreases exponentially in (s — T c )/e, which 
implies that the area between them is at most of order so* 2 ' 3 . Thus ( |5.21 ) holds for 
H > 0(ea 2 ' 3 ). But for smaller H, it is trivially satisfied. 

We can summarize the properties obtained so far in the following way. 

Proposition 5.6. There exist constants C , n, h\, hi, hs > such that 

T 4 /3 /(e | lo g CT |) 



F t0 ' X0 {A -A< -H) ^ - e - KH3/2 / a2 +- e _wr " 

c — ■ 

F to ' x °{A-A> +H} ^ -e- KH ^ a2 ^^+-l [0Ma4/3) (H)e- KH/{£llogal) 
for ^ H ^ /iic 2 ' 3 |log(r|. In addition, for all H ^ /13 we have 



(5.22) 
(5.23) 

(5.24) 



As an immediate consequence, we obtain the following estimates on the moments of 
the deviation of the area. 

Corollary 5.7. There exist positive constants C, c\, c<i such that 

E t0 ' X0 {A(e,a) -A} ^ -Ca 4 / 3 |loge| 2 / 3 (5.25) 

E t0 ' X0 {A(s,o-)-A} ^C[(eVa 2 |loge|)|logCT|] (5.26) 

E t0 ' X0 {(A(e,a) -A) 2 } ^ C((j 4/3 |loge| 2/3 ) 2 (5.27) 

provided c\E < <7 4 / 3 /|logcr| 2 and cr 2 ' 3 |logcr| ^ C2/|loge|. 
Proof: By partial integration, 

/■OO /'OO 

E t ,x 1_4 _ _4| = / F to,x {A-AyHJdH- / P 4o ' x '° {A - A < -H) dH. (5.28) 
Jo Jo 

The first integral be evaluated by splitting it at H = C(e V <r 2 |loge|)|logcr|, /i2cr 4 / 3 and 
/13, and the second one at Ca 4 ' 3 1 log e| 2 ' 3 and /13. Each time, the integral over the first 
interval dominates. The estimate ( 5.27[) is obtained similarly. □ 
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6 The large amplitude case 



We consider finally the large amplitude case ao ^ 71 £, but with noise intensity a satisfying 
a 2 < ^/oq£. By symmetry, we may again concentrate on a half-period [—1/2,0]. Let xt 
be the solution of the SDE ( 2.1 1] ) starting at time to = —1/2 in the upper well, i. e., near 
x* (t). Recall that the solution xf ct of the deterministic equation ( |3,1| ) with the same 
initial condition tracks x+(t) until time —t c , and jumps to the other potential well at 
x*_(t) after a delay of order £ 2 ' 3 a . 

We introduce a time t* = —t c + cq£ 2 ' 3 a just before the jump. Then we know that 
for to < t = —t c + s < £*, 



det _ 



1 



at) 



a(t) x \s\ V aj 74 ^ V (e^) 1/3 , 



(6.1) 



compare ( |3.18| ) and (|4.6|) . The fact that xf et behaves in this way follows from the fact 
that x%(t) - x c x a/M^) dominates x*{t) + - xf et for t ^ -t c - e 2 ^^ 1/6 , c.f. (|I|). For 
larger t, we know from ( |3.16 ) and ( 3.1 7| ) that xf et — x c x (e^/ao) 1 ' 3 . 



6.1 The transition time 

Let us again start by investigating the distribution of the stopping time 

t° = inf{s > to- x s = 0}. 
The following result shows that r° is likely to be close to t*. 

Proposition 6.1. There exist constants C, c\, k > such that 
• for t ^t < t*, 



P *0,*0{ T << } <C r(Lj0 + J ] ( . xp 



(6.2) 



a 2 CW 2 



(6.3) 



/or t* ^ t ^ f* + ciy/a~o, ^ #t* ^ p and any p > x c 



vt*,x t + / 



{r° >t} < 3exp 



£y/a 



a- 



% ,b (t-t* 



|log(e 2 a )| Vlogp e 2 / 3 



1 



(6.4) 



Proof: First note that ( |6.3[) is a direct consequence of Proposition L2 as we are in the 
stable case. 

In order to prove ( |6.4[) , we consider again the stochastic process yt = xt — x c , satisfying 
the SDE 



dy t = - [p(t + t c ) - VSy 2 - y 3 ] dt + 4= dW t , 

e ye 



y t * = z t * - x c 



(6.5) 



where ( |3.14D implies p(t + i c ) ^ — e 2 ' 3 a for £* ^ £ ^ t* + ciy^. Note furthermore that 
— V3y 2 — yf ^J — y 2 for i ^ r°. By Gronwall's inequality, it follows that yt ^ 2^ for i ^ r°, 
where zt is defined as the solution of the time-homogeneous SDE 



£ v £ 



z t * = x t * - x c 



(6.6) 
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For any 5q > 0, we can write 



{r° > t] < P'*-***/ sup z s > (J ) + P* Vt *{-x c < z s ^ 5 Vs G [t*,i]}. (6.7) 



Since ( p,6|) is an autonomous SDE, it is easy to see that the first term on the right-hand 
side can be bounded by 



r>t* ,x 






{ sup z s > 5 \ ^ CT-— + l) e"^ 2 , (6.8) 



which can be made as small as we like by taking 5q sufficiently large. In order to estimate 
the second term, we introduce A = ce 2 ' 3 a , where c > 1 will be chosen later, and define 

Q = sup P°' 2 ° { -x c ^z s ^5 Vs£ [0, A] } . (6.9) 

Using time homogeneity and the Markov property, we can write 

p**>*t*{_a. c ^ Zs <; s Vs G [t*,t]} ^ g(*-**)/A-i. ( 6-10 ) 

The result is thus proved if we manage to bound Q by a term exponentially small in 
£^/a ~/a 2 . 

In order to estimate Q, it is convenient to introduce the process z\ = —zt/(e 1 ' s a ), 
which obeys the SDE 

ds^Ii + fla-.^dw,, * = -^- («D 

Let /9 = x c /(e ' s a ) and S = 5o/(e 1 ' 3 a )■ Using again Markov property and time- 
homogeneity shows that Q ^ Qi + Q2 + Q3, where 

Q 1 = sup F°'~ zo {z s < -1 Vs G [0, A/3]} 

-<5<zo<-l 

Q 2 = sup F°' z °{z s < 1 Vs G [0,A/3]} (6.12) 

Q 3 = sup P°' 5 ° { z s < p Vs G [0, A/3] } . 

Since 1 + z 2 ^ 1 V|£:|, each term can be easily estimated by comparison with an appropriate 
linear or ^-independent equation. Consider for instance Q\. We know that zt lies above 
the solution 2® of the linear SDE 

1/6 

di *° = ~%~ z " dt ~ ^fr/6 dw <> $ = ~ 6 > ( 6 - 13 ) 

x 



£ 2/3 * ^/e 1/6 



the solution of which at time A/3 is a Gaussian random variable with mean — 5e c ' 3 and 
variance (1 — e _2c ' 3 )<7 2 /(2ey/ao). We can thus estimate 

<h < P».-{li /3 < -1} < e,p{-£^il^£}, (6.14) 

provided 5e _c ' 3 < 1, i.e., c > 3 log 5. Now, Q2 and Q3 allow for similar bounds, and the 
result thus follows from Q6.8| ) and ( |6.10|) , taking 6q and c sufficiently large. □ 
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6.2 The case Ha 



We now examine the process Yt defined in (4.20), which describes deviations from the 
deterministic area. Using Lemma |4.11| and (3.1), it is easy to check that 



ri(*Vo) 



lloefe 



2/3 -lM 



and r(i*,i )-l, 



(6.15) 



where to = — 1/2- Applying Remark \LQ and Corollary 4.10) , it is straightforward to check 
that the distribution of Yt is close to a Gaussian with variance proportional to a 2 e. 

The situation changes, however, for t > t*, because the deterministic solution crosses a 
zone of instability between x c and — x c , compare (|Q|). This instability causes a spreading 
of paths which we will now analyse in more detail. Let us introduce times t\ and t\ such 
that 



„dct 






ci(e^/o^) 



1/3 



„dct 



S, 



(6.16) 



where c\ > and 5 < x c . Then t* — t* and t\ — t* are both of order e 2 ' 3 a . We 
now proceed to determining the behaviour of C(*)> defined in (4.5), which measures the 
spreading of paths around the deterministic solution. 



Proposition 6.2. Let zt = 

C(t) X (ey^) 

C(t) x (eV5o) 



„det 



Then there exist constants C, K > such that 



-1/3 



-5/3 A 



C(<) < C[(eVo?) 



-5/3 -K(t-t* 2 )/e 



+ 1] 



for t* < i ^ i* 


(6.17) 


for ti^t^tl 


(6.18) 


for t% < t ^ 0. 


(6.19) 



Proof: ( 6.17| ) follows by an elementary calculation. Next, consider the time interval 
t* ^ t ^ #2- ^he var i a ble z^ satisfies the differential equation 



¥ = ;W« + ^ 



"t > 



fj,{t) = A c - Acos(2irt). 



(6.20) 



Note that ^ is monotonously decreasing and n(t) ^ for the times under consideration. 

2v3zt — 3z?. It follows that for 



The linearization of the drift term F at zt is a(t) 

zt 2V3z - 3z 2 



q^s^t^q, 



a(t, s) 
and thus 



a(u) du ^ e 



V3z 2 



dz = e log 



V3- 



zt 



z%(y/3- z s ) 



e a{t,s)/e ^ ft 



(V3- 



ZfJ 



>2' 



2;f(\/3-z s ) 
More careful estimates, based on the inequalities 

show that e a ( t,s >' £ is also bounded below by a constant times (zt/z s ) 2 . 
be computed in the same way, by performing the change of variables s 



(6.21) 



(6.22) 



(6.23) 

Now £(£) can 
-* z s , yielding 
(6.18). Finally, ( |6.19 ) follows easily from the fact that we are again in the stable case for 



t ^ q. 



D 
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We now return to the SDE (gj) for a(t) given by ( gig) and 6(y,i) = F(y,i) - 
a(t)y. Following the proof of | BG2| , Proposition 3.10], it is easy to establish ( |4.9Q for 
all /i ^ h()(£^/ao) 5 ' e and all t. The condition on h stems from the fact that the linear 
term a(t)yt should dominate the nonlinear term b{yt,t) for all realizations to satisfying 

The condition on h implies that we need to require a ^ (ey^) 5 ' 6 for (|4.9| ) to be 
of interest. Then the maximal spreading of paths will typically be of order (JyC^I) x 

-5/6 



o-(£Vao 



Since |a'(t)| is no longer bounded for t* ^ t ^ ££> we cannot apply Lemma 4.11 to 
compute the integrals ( 4.25 ) and ( (4.26 ). However, using the same change of variables as 
in the proof of Proposition |6.2j, it is not difficult to establish that 



(6.24) 



ri(0,t )^(ev^) _2/3 and r(0,t )^e _2/3 ao /6 , 



where again to = — 1/2. 



The following proposition now follows immediately from Corollary 4.10 



Proposition 6.3. Assume that a ^ (e-^/ao) 5 ' 6 . There exists a constant hi such that for 



^{\Y \ >H}^ -exp{-K ^_ u/3 }- 



(6.25) 



Finally, Proposition 47 can also be applied to show that for cr|loge| ^ const (et/oq) 5 '®, 

a 2 1 log e\ 



E 



to,a:o 



M|=0 



(£v^) 2/3 
Var(y )^^ 2 (ev^) 1/3 - 



(6.26) 
(6.27) 



6.3 The case lib 

For a larger than (e-^/aTj) 5 ' 6 , the strong dispersion of trajectories near time t\ prevents us 
from applying methods of Section |j. However, methods similar to those of Section || can 
be applied to obtain some information. 

The following result allows to estimate deviations of Yj* in a larger domain than Propo- 



sition 4.9 



Proposition 6.4. There exist constants C , k, h± > such that 



C 
F to ^ iY t * < -H\ ^ -e 

1 ' e 



-k(H 3 / 2 \/e^)/* 2 



F to > xo {Y t * >+H} <e 



-nH 2 /(a 2 e) 



C 

+— e" 

e 



-Ktf/(<7 2 ri(t*,t )) 



(6.28) 
(6.29) 



forOs^H < /ii(e y / ?io") 1 / 3 ri(t*,to), where ri(t*,*o) - |log(e 



2/3„-V6^ 



Proof: The proof is almost the same as the proof of Proposition 5.4, the only difference 
lying in the different behaviour of C(*)> gi ven i n Q6.1[ ), which requires to distinguish between 
T + t c ^ — i/ao, — y / oo ^ t + i c Sj — e 2 ' 3 a , and the remaining r up to t*. D 
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Proceeding as in the proof of Proposition 5.5 , but using Proposition 6T for the tran- 
sition time, we obtain 

Proposition 6.5. There exist constants C, k, hi, h% > such that for all xt % £ [—L,L] 

(L x I), and h 2 (e^) 2 /3\\og(e^a- )\ ^H^h 3 , 

p*Vt* |y < _#} ^ ^Z e -*H/a* + e -KH^/{^e) (6.30) 

pi*,a; f * jy > _|_g-"l ^ _ e -K(e v /s o) 1/3 - H '/( " 2 l lo g( £ v /s o)l) . (6.31) 

Corollary 6.6. For h 2 {e^) 2 / 3 \log(e^)\ ^ h^e^) 1 ^^, to), 

F t0 > X0 {Y < -H} ^ — e~ KH3/2 /^ (6.32) 



pio.zo jy Q > +JJ\ ^ ^ e -K(e v ^5') 1/3 ^/((T 2 |log(e v ^')|) _ (6.33) 



The required lower bound on H only allows us to conclude that expectation and 
standard deviation of Yq are smaller than a constant times (ey / ao) 2 ' 3 |log(ey / ao)|, although 
the above estimates are already very small for H = /i2(ey / ao) 2 ' 3 |log(ey / ao)|. 
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